Bose-Einstein Quantum Phase Transition in an Optical Lattice Model 
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Bose-Einstein condensation (BEC) in cold gases can be turned on and off by an external potential, 
such as that presented by an optical lattice. We present a model of this phenomenon which we are 
able to analyze rigorously. The system is a hard core lattice gas at half-filling and the optical lattice 
is modeled by a periodic potential of strength A. For small A and temperature, BEC is proved to 
occur, while at large A or temperature there is no BEC. At large A the low-temperature states are in 
a Mott insulator phase with a characteristic gap that is absent in the BEC phase. The interparticle 
interaction is essential for this transition, which occurs even in the ground state. Surprisingly, the 
condensation is always into the p = mode in this model, although the density itself has the 
periodicity of the imposed potential. 



I. INTRODUCTION 

One of the most remarkable recent developments in 
the study of ultracold Bose gases is the observation 
of a reversible transition from a Bose-Einstein conden- 
sate to a state composed of localized atoms as the 
strength of a periodic, optical trapping potential is var- 
ied 12 1|. This is an example of a quantum phase tran- 
sition |3( where quantum fluctuations and correlations 
rather than energy-entropy competition is the driving 
force and its theoretical understanding is quite challeng- 
ing. The model usually considered for describing this 
phenomenon is the Bose-Hubbard model and the transi- 
tion is interpreted as the superfluid-insulator transition 
that was studied in £| with an application to He 4 in 
porous media in mind. The possibility of applying this 
scheme to gases of alkali atoms in optical traps was first 
realized in Q • The article [(| reviews these developments 
and many recent papers, e.g., 0, H, 0, E3, EE E3, 113, Q 
are devoted to this topic. These papers contain also fur- 
ther references to earlier work along these lines. 

The investigations of the phase transition in the Bose- 
Hubbard model are mostly based on variational or nu- 
merical methods and the signal of the phase transition is 
usually taken to be that an ansatz with a sharp particle 
number at each lattice site leads to a lower energy than 
a delocalized Bogoliubov state. On the other hand, there 
exists no rigorous proof, so far, that the true ground state 
of the model has off-diagonal long range order at one end 
of the parameter regime that disappears at the other end. 
In the present paper we study a slightly different model 
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where just this phenomenon can be rigorously proved and 
which, at the same time, captures the salient features of 
the experimental situation. 

The model is that of a hard core gas on a cubic lat- 
tice at half-filling (i.e., when the particle number is half 
the number of sites). The 'optical lattice' is modeled 
simply by a periodic, one-body potential A(— l) x , where 
(— l) x = +1 on the A-sublattice and (— l) x = — 1 on the 
B-sublattice. Thus, the Hamiltonian, expressed through 
bosonic creation and annihilation operators, equals 

H = ~lY.^ a y + a x«y)+ A E(- 1 ) Xa x a x 

<xy> x 
+C / 5I a x a x(4ax ~ !)• (!) 

X 

The sites x are in a d-dimensional hypercubic lattice, 
and (xy) stands for pairs of nearest neighbors. The case 
considered in this paper is the hard-core interaction U = 
oo. Apart from the periodic potential, this is also the 
Bose-Hubbard model with infinite on-site repulsion. 

Note that in our model the troughs of the optical po- 
tential correspond to the B-sublattice where the periodic 
potential is negative, and the crests correspond to the A- 
sublattice. Often in the Bose-Hubbard model the whole 
lattice itself is used to approximate the troughs alone. 
Roughly speaking, half-filling in our model corresponds 
to filling factor 1 in the Bose-Hubbard approximation. 

As is well known, the model with U = oo can also 
be viewed as the XY model of a spin 1/2 system |l5). 
The periodic potential then corresponds to a staggered 
magnetic field. This will be explained in the next section. 

We are able to prove the following facts rigorously for 
3 or more dimensions. The A-T-phase diagram at half- 
filling (e.g., mean density g = ^) is shown schematically 
in Figure ^ 
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1. If both A and the temperature T are small, then 
there is Bose-Einstein condensation (BEC). In this 
parameter regime the one-body density matrix has 
exactly one large eigenvalue (in the thermody- 
namic limit), and the corresponding condensate 
wave function is y(x) = const. 

2. If either T or A is big enough, then the correla- 
tion function (the one-body density matrix) decays 
exponentially, and hence there is no BEC. In par- 
ticular, this applies to the ground state (T = 0) for 
A big enough. 

3. The Mott insulator phase is characterized by a gap, 
i.e., a jump in the chemical potential (at zero tem- 
perature). We are able to prove this, at half-filling, 
for big enough A. More precisely, there is a cusp in 
the dependence of the ground state energy on the 
number of particles; adding or removing one par- 
ticle costs a non-zero amount of energy. We also 
show that there is no such gap whenever there is 
BEC. 
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FIG. 1: Schematic phase diagram at half-filling 

4. The interparticle interaction is essential for items [3 
and|3 Non-interacting bosons always display BEC 
for low, but positive T (which depends on A, of 
course) . 

5. For all T > and all A > the diagonal part 
of the one-body density matrix (a x a x ) is not con- 
stant. Its value on the A-sublattice is constant, 
but strictly less than its constant value on the B- 
sublattice (for a finite system with periodic bound- 
ary conditions) and this discrepancy survives in the 
thermodynamic limit. In contrast, in the regime 
mentioned in item ^ the off-diagonal long-range 
order is constant, i.e., (a x a y ) ~ const, for large 

|x-y|. 

We give explicit expressions for the curves, sketched 
in Figure ^ defining the regimes for which the above 
statements are proved here (see Eqs. © and <|13[) below.) 



We are not able to make a rigorous statement about the 
intermediate regime, but we believe that there is only a 
critical line separating the BEC and the Mott insulator 
phases. 

We focus here on lattice dimensions d > 3 but, us- 
ing the technique employed in |l6j . an extension to the 
ground state in two dimensions is possible. Other possi- 
ble extensions are mentioned in the next section. 



II. DETAILED DESCRIPTION OF THE MODEL 

We write the Hamiltonian (Q of the lattice gas with 
U = oo in terms of the creation and annihilation opera- 
tors, a x and a x , for particles at lattice site x £ A, with 
A a finite hypercubic lattice with L d sites, L being an 
even integer. We impose periodic boundary conditions. 
Because of the hard-core condition, there is at most one 
particle at each site, and thus the creation and annihila- 
tion operators can be represented as 2 x 2 matrices with 
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for each x G A. The correspondence with the spin ma- 
trices 

*-K!i)-*-K!o)-*-Ki- i) 

is 

a x = i§ x + iS x = S£ , a x = S x — iS x = S x , 

and hence a x a x = S x + \ . Adding a convenient constant 
to make the periodic potential positive, the Hamiltonian 
for U = oo is thus equivalent to 

h = -iJ2(s+s; + s x s+) + \J2 [h + (-i) x ^x] 

<xy) x 

= -E(^y + ^ y ) + A E[l + (- 1 ) XS 'x]-( 2 ) 

<xy> x 

As explained in the introduction, (— l) x = ±1 on alter- 
nating sites. Without loss of generality we may assume 
A > 0. Note that the subtraction of the 'diagonal' terms 
in the kinetic energy has the effect of a chemical poten- 
tial and as a consequence the unique ground state of PJl 
has particle number N = s|A|. We postpone the proof 
of this assertion to Appendix lAl 

The presence or absence of Bose-Einstein condensation 
is expressed through the reduced one-particle density ma- 
trix 

7(x,y) = (a x a y ) 

where (•) denotes the expectation value in the thermal 
equilibrium state or the ground state considered. BEC 
occurs (by definition) if the L d x L d matrix 7(x, y) has 
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an eigenvalue of order N in the thermodynamic limit 
A — > oo, N — » oo, with g = N/\A\ fixed. 

We shall prove that for g = 1/2 and d > 3 the thermal 
equilibrium state of J2J) shows Bose-Einstein condensa- 
tion for small A and low temperature T, while for large 
A or T the condensation disappears. For d — 2 this is 
true only in the ground state. Here, g stands for the 
average density, since we are using the grand-canonical 
ensemble where the particle number is not fixed. Note 
that we are not dealing here with a dilute system and 
the condensation is always depleted, even in the ground 
state. 

We remark that the Hamiltonian is invariant un- 
der the following two unitary transformations, which will 
be used throughout the paper. 1 ) Uniform rotation 
around the S^-axis, in particular the map — > — and 
Sx — * —S% at all sites. The corresponding conserved gen- 
erator is the total particle number. 2) The particle-hole 
symmetry, which corresponds to: «-» ~* ~ 

at all sites, followed by a unit-vector translation in any 
of the lattice directions. For the latter symmetry the 
half-filling is essential. 

Our analysis of the system proceeds via the following 
steps. In the next Section lTTTI we prove infrared bounds on 
the two-point function (S^S^p + S^S^p) in momentum- 
space. The essential ingredient here is reflection positivity 
of the Hamiltonian J3J) and the closely related property of 
Gaussian domination [Tt| . The bound obtained depends 
on A and the temperature. A sum rule fixes the sum over 
p of the two-point function and the infrared bound leads 
to the conclusion that for small A and low temperatures 
the contribution from p = remains non-vanishing in 
the thermodynamic limit. This proves BEC. 

Reflection positivity is essential for our proof of BEC 
and forces the periodic potential to have period 2. Other 
generalizations can be accommodated, however, such as 
a more general lattice in which we add hopping to next 
nearest neighbors or the addition of nearest neighbor 
interparticle repulsion |17| . For simplicity we concentrate 
here on the simple cubic lattice with on-site repulsion 
only. 

In Section IIVI we show that the existence of BEC im- 
plies the absence of an energy gap for adding or removing 
a particle, and that the energy viewed as a function of the 
density has a unique tangent at Q = \, i.e., the chemical 
potential is continuous. 

The absence of BEC for large A or high temperatures 
T is proved in Sect ion E] The technique applied here is a 
path space representation of the two-point function that 
follows from the Trotter product formula. This represen- 
tation allows us to derive exponential decay of the two- 
point function, provided A or T are sufficiently large, and 
hence absence of long range order. The magnitude of A 
or T enters through the suppression in the path space in- 
tegral of long contours connecting a pair of lattice points 
if these parameters are large. The same representation is 
used for proving the existence of a gap, as explained in 
item in the introduction. This method is quite robust 



and easily extends to a periodicity of the optical lattice 
potential different from 2, for instance. 

In Section IVII we show that the particle density (or 
3-component of the spin) oscillates with the period of 
the staggered field if A ^ 0, in contrast to the conden- 
sate wave function which is independent of A and x. For 
this effect the interaction is essential, as remarked in Sec- 
tion Em Without the interaction there is always BEC 
(for low T) and the condensate wave function is never 
constant (for A ^ 0). 

In Appendix 1X1 we shall prove that the ground state 
of J2)l has total spin 0, which in the lattice gas language 
means that the lowest energy of @ is obtained when 
the particle number is |A|/2. The essential ingredient in 
our proof is again reflection positivity of the Hamiltonian 
. We also show that the canonical partition function 
is maximal at half-filling. 



III. PROOF OF BEC FOR SMALL A AND T 

In this section we are going to show the occurrence of 
BEC for small A and low enough temperature. The main 
result is the following. 

THEOREM 1 (Existence of BEC). Let E p = 

Yli=iO- ~ C0S (Pi)) (where pi denotes the components of 



p) and 



1 



Cd = 



(2ir) d 

In the thermodynamic limit 



&jXp £ 7(*,y) 



-fir. 



1 1 



x,yeA 



1/2 



1/2 1 



with f3 = l/(k-QT) the inverse temperature. More- 
over, if ip(x.) denotes the (normalized) eigenfunction 
corresponding to the largest eigenvalue o/ 7(x, y), then 
liniA^oo |A| _1 | ^ x <f?(x)| 2 = 1, implying that the con- 
densate wave function is constant in the thermodynamic 
limit. 

Note that Cd is finite for d > 3. Since the largest 
eigenvalue of 7(x, y) exceeds | TV | 1 J2 X y t(x, y), BEC is 
proved if the right side of the expression is positive. 
This is in particular the case, for large enough /?, as long 
as 

a 1 d(d+l) 

In d = 3, C3 » 0.505 01 > an d hence there is BEC for 
A < 0.960. In [13 it was also shown that ded is monotone 
decreasing in d, which implies a similar result for all d > 
3. 
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The main tool in our proof of Theorem^s an infrared 
bound as in ^tJ- The statement is as follows. For A and 
B bounded linear operators, denote by 



(A,B) = J Tr (Ae~' pB Be-^-'^ a ) ds/Tre 



0H 



the Duhamel two-point function. For p € A* (the dual 
lattice), p ^ and Sp = | A!" 1 / 2 J2 X Si exp(ip • x), we 
claim that 



(sp,si p )< 



20E P 



(4) 



The same is true with Sp replaced by S 2 . This inequality 
will allow us to prove the bound © above. Moreover, 
we will infer from @ the fact that there is only one 
large eigenvalue (of order |A|) of the one-particle density 
matrix, and the corresponding eigenfunction is constant 
(in the thermodynamic limit). 

We start by proving Q. The main ingredient is Gaus- 
sian domination. More precisely, let h be any real-valued 
function on A, and 

Z(h) = Tr cxp \-0K(h)} , 

where K(h) is the modified Hamiltonian 

K (h) = £(|(£i-4-fc* + /iy) a -.s2s£) 

(xy) 

+A^(-1) X ^. 



Note that for h 
constant. 



this operator agrees with H up to a 



Lemma 1. For all real-valued functions h, 
Z{h) < Z(0). 

Proof. We perform a unitary transformation that takes 
S 2 i-> — S 2 and S 3 h-> — S 3 on the B-sublattice. Since the 
trace does not change under unitary transformations, we 
have Z(h) = Tr exp[-(3K(h)} with 



K(h) 



Ev, 

(xy) 



l(Si -S y -h x + h y ) 2 + S^)|AVS 



Compared to K(h), the sign in front of the S 2 S 2 term 
has changed, and the (— l) x has vanished. The operator 
K(h) thus obtained is translation invariant. Since S 1 and 
S 3 are real self-adjoint matrices, and S 2 is imaginary and 
self- adjoint, we meet exactly the conditions for applying 
the result in [tH Lemma 6.1 and proof of Theorem 4.2] 
to prove the lemma. □ 

The infrared bound J3J follows from this lemma by 
using the negativity of the second derivative, 



< 0. 



By performing the derivative, we obtain 



(xy) 



where A = E( X y>( 5 'i ~ S y)( h x ~ h y)- We proved this 
inequality only for real- valued h, in which case A = A\ 
but it automatically extends in a standard way |17| to 
complex-valued h. In this case note that the adjoint 
of A agrees with its complex conjugate. Now, choos- 
ing /i x = exp(ip ■ x), we obtain By invariance of 
the Hamiltonian under rotations around the S 3 -axis, the 
statement is also true with Sp replaced by S 2 . 

We now want to use 1 7. Theorem 3.1] to relate the 
Duhamel two-point function to the ordinary thermal two- 
point function. For that purpose, we have to evaluate the 
double commutators 

[5i,[i/,5i p ]] + [5 2 ,[i/,5 2 p ]] 

= -m( H - ^ A|A| + 2 E s * s y cosp • (x - : 



(xy) 



Let C p denote the expectation value of this last expres- 



C p = ([Sp, [H, Sip}} + [S 2 , [H, S 2 _ p }}) > 0. 

The positivity of C p can be seen from an eigenfunction- 
expansion of the trace. From 0, Corollary 3.2 and The- 
orem 3.2] and @ we infer that 



s 2 s 2 




£ coth J(3 2 CpE p /4. (5) 



Using cothx < 1 + 1/x and Schwarz's inequality, we ob- 
tain for the sum over all p ^ 0, 



< 



-V— ifV — 

E„ + 2 V 2-. Er, 



1/2 



p^O 



p#0 



p#0 



1/2 



We have J2 P eA> C p = - 2 ( H ) + MMi which can be 
bounded from above using the following lemma. Its proof 
follows exactly the same lines as |Tj], Theorem C.l]. 



Lemma 2. The lowest eigenvalue of 

2d ^ 2d 

■ 2 xE 5 y* _ 2'' 



^ 2d ^ 2d 

' o S * E ^y* - o^x E S yt + AS * 



(6) 



is given by -i[d(d + 1) +4A 2 ] 1/2 . 



Since the Hamiltonian can be written as a sum of terms 
like 10, with y-j the nearest neighbors of x, we get from 
this lemma the lower bound 



H > 



|A| 



[d(d + 1) + 4A 



21 V2 



±A|A|. 
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With the aid of the sum rule 
peA* 

(which follows from (S 1 ) 2 = (S 2 ) 2 = 1/4), we obtain a 
lower bound, in the thermodynamic limit, 

Mm ^-AS^ + S'S 2 ) 

A-*oo |A| 

>l-l(kAd(d+l)+^ 2 ] 1/2 c d ) 1/2 -^c d . (7) 

The connection with Bose-Einstein condensation is as 
follows. Since H is real, also 7(x, y) is real and we have 

T (x,y) = {S x + 5 y > = (^+Sx^ 

Hence, if tpo = | TV | 1 / 2 denotes the constant function, 

(< Po \ 7 \<p } = (S 1 S 1 + S 2 S 2 ), 

and thus the bound Q) implies (J3J. In addition we claim 
that the infrared bounds imply 

(<p|7|<P> < const. |A| 2/d 

for any normalized (p that is orthogonal to ip^ 1 with a 
constant that is independent of cp. To see this, consider 
the positive definite matrix (S+S" }, with p ^ 0, q ^ 
0. The infrared bound (JSJ implies that the diagonal of 
this matrix is bounded by (SpSZ p ) < const. |p|~ 2 < 
const. |A| 2 / d for p / 0. Moreover, the matrix is almost 

diagonal in the sense that (Sp~ SC q ) ^ on ly if Qi — Pi or 
qi = Pi in (by invariance under translation by two lattice 
sites). The largest eigenvalue of such a matrix is bounded 
above by 2 d times the maximum on the diagonal, namely 
2 d const. |A| 2 / d < |A|. This proves our claim. 

We conclude that 7(x, y) has exactly one large eigen- 
value, with corresponding eigenfunction equal to ipo as 
A| — > oo. I.e., the condensate wave function is constant. 
This is in contrast to the particle density, which shows 
the staggering of the periodic potential. We show this in 
Section IVll below. It also contrasts with the situation for 
zero interparticle interaction, as discussed in Section lYlII 

IV. NO CUSP, NO GAP 

The system is in a Mott insulator state at zero tem- 
perature if a finite change in the chemical potential is 
required to change the particle number in the ground 
state. We refer to this as a gap in the chemical poten- 
tial. More precisely, if Ek denotes the lowest energy of 
@ restricted to the sector of i|A| + fc particles (which 
corresponds to S^ ot = ^ x S%. = k in the spin language), 
a gap means that, for all fc, 

E- k + E k - 2E > c\k\ (8) 



for some c > independent of A and k. Note that 
particle- hole symmetry implies that E-k + Ek — 2Eq = 
2(E W -E ). 

In the next section we prove © for sufficiently large A. 
In this section we will show that whenever there is BEC 
then JSJ fails. In fact, we will prove that 

E k - E < (9) 

for some fc-dependent Ck > 0, which is independent of A, 
i.e., that (JSJ does not hold for finite k. 

This, however, does not rule out the possibility that 
the macroscopic system still acts as an insulator. To 
show that this is, indeed, also not the case, we prove 
that fails for macroscopic k as well. More precisely, 
we will show that the thermodynamic limit of the ground 
state energy per site, eoo(o) = limA— «x> Ek/\A\, where 
k = (g — |)|A|, satisfies a bound 

< e^g) -eood) < const, (g - \) 2 (10) 

for g close to 1/2, i.e., that there is no macroscopic cusp 
in the energy at half-filling. 

We will first prove (|1(J|) . With |0) being the ground 
state of H , and with y some point in the lattice, consider 
the states 

h/VHe- s -(4 + i)|0>. 

The motivation is the following: we take the ground state 
and first project onto a given direction of 5* 1 on some 
site y. If there is long-range order, this should imply 
that essentially all the spins point in this direction now. 
Then we rotate slightly around the S* 2 -axis. The particle 
number should then go up by e|A|, but the energy only 
bye 2 |A|. 

The norm of |^ y ) is given by 

<Vvh/v) = (o|s y + i|o> = i, 

where we used the symmetry S 1,2 — > — S* 1,2 . We want 
to find an upper bound to the average energy of these 
states, more precisely, 

' ' y 

Here Eq denotes the ground state energy of H, which is 
obtained at half-filling (see Appendix Ia"|> . We claim that 
the inequality 

e -i^ t 2 otiJe ieS t 2 ot < H + i£ [ Hj sf ot \ + const. e 2 |A| (11) 

holds for some constant depending only on d and A. To 
see this, consider, for self-adjoint operators A and C, 

F A (e) = e ieC Ae- ieC . 

By Taylor's formula 

F A (e) < F A (0) + eF^O) + \e 2 sup \\F^ V )\\. (12) 

0<r]<e 
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Note that the last norm is actually independent of r], Since there is BEC, ({S tot ) k (S^ ot ) k ) > c' k \A\ 2k for some 
since e~ leC is unitary, and is given by the norm of the c' k > modulo lower order terms as |A| — > oo. All we 

have to show is that [(S^ ot ) k , [H, (5+ t ) fc ]] < 4\A\ 2k - 1 for 
some constant c k . This is clear, however, since altogether 
there are |A| 2fc+1 factors, and the 2 commutators reduce 
the power by two. Hence we obtain JjJJ. 



double commutator [C, [C, A]]. After evaluating the dou- 
ble commutator for the case in question, a simple bound 
gives l(TT1) . 

Consider now the first term on the right side of 
We obtain 



y 



(0\( Sy + ±)(H-E )(S y + ±)\0) 
2SlS 2 y -\J2(-l) y S 3 y 



<xy) 

= -§(£fc-§A|A|), 

where we used rotational symmetry in the last step. The 
second term, 



is zero by symmetry, as can be seen in the following way. 
The diagonal terms are zero by the symmetry S^' 2 — > 
— 5* ,a at all sites. The off-diagonal terms are 

( |£tot[- ff >'Stot] + [-ffj'S'totl'S'totl )! 

which is zero by the symmetry S^ 3 — > — S^' 3 at all sites, 
followed by a unit- vector translation in any of the lattice 
directions. We therefore get that 



AE < const. e 2 |A| 



|£o| + ^A|A| 
|A| ■ 



It remains to evaluate the average particle number of 
the states considered. Using that Sf ot |0) = 0, we obtain 



_^(0|^ + i)e- 



^sine<0|(^ ot ) 2 |0), 



which is of the order e|A| if there is BEC. Choosing e 
proportional to the chemical potential /i, we obtain as an 
upper bound for the ground state energy of H — /iSf ot 



£>o(l-|A|- 



CM 2 |A| 



for small fi, with c > in the case of BEC. By taking a 
Legendre transform, we arrive at (|10(l . 

To prove © we use as a trial state (S' t ^ t ) fc |0), with 
k > 1. Using the particle- hole symmetry of H as well as 
the fact that we are considering the ground state, we get 
the bound 

P <T? , i ([(sr ot )Mff,(s+ t ) fc ]]) 

Ek£Eo + 2 ((s^ns+ t ) k ) ■ 



V. GAP FOR LARGE A, AND ABSENCE OF 
BEC FOR LARGE A OR HIGH TEMPERATURE 

We shall now present explicit bounds for a region of 
values of (A, T) for which BEC is absent. This region 
includes: 

(i) all A > at k B T > d/(21n2), 

(ii) all T > at A > such that A + |e(A)| > d. 

The absolute (i.e., without specifying N) ground state 
energy per site for a finite A, which is always obtained at 
half-filling, is denoted by e(A). Note that e(A) < 0. 

In this regime, the particles are localized in the 
sense that the transition amplitudes decay exponentially. 
Short excursions occur locally in "space-time" , however 
a long distance transition requires a linked chain, or per- 
colation, of such local events and the amplitude for that 
decays exponentially, as in sub-critical percolation mod- 
els. One may discern here two distinct mechanisms con- 
tributing to the localization: at high A localization is 
caused by the confinement due to the staggered struc- 
ture of the potential, whereas at high temperatures it 
is a combined effect of the exclusion (no more than one 
particle at a lattice site) with the reduced amplitude for 
coordinated moves by neighboring particle. The above 
picture is made precise in a representation of the matrix 
elements of e~^ H which in effect involves an "imaginary 
time" . 

A more inclusive statement of the condition under 
which our results hold is 



(13) 



where / = /(/3, A) = — (/5| — 1 lnTre"^ is the free 
energy per site, which satisfies 

-/(/3,A) >max{|e(A)|,/3- 1 ln2-iA}. 

Under condition l|T3|l. we define v > by 

d 



A-/ 



(l-e-K x -fi) < 1. 



THEOREM 2 (Mott insulator phase). Throughout 
the regime where ijj^jj holds the thermal average two-point 
function decays exponentially. More specifically, for any 

Tr ata„e~@ H , 
7(*,y)= T ; e % H <Cfee-gl-*l (14) 
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with C(: — [1 — "] 1 . Similar decay also holds in 
the finite volume ground state ( corresponding to the limit 
/?->• oo in gpj. 

Moreover, the ground state energy Ek for particle num- 
ber i|A| + k satisfies 



E k + E_k - 2E Q = 2(E k -Eq) > c \k\ 



(15) 



with c = 2(A + |e(A)| — d), which is strictly positive for 
large A (independently of the volume \K\). 

Eq. I|15|l is, in fact, the consequence of a more explicit 
result: If Vk denotes the projection onto the subspace of 
fixed particle number N = 4|A| + k, then 



TrV k e- 0H 
TrV e-P H 



< 2 e -«/3|fe| 



|A| 



n \k\ 



B{a) \k\ 



(16) 



for any a > for which 







B(a)=d / e -^+l/o(ftA)l-°)* t ft<i j 
J o 

with /o(/3, A) denoting the free energy for the system with 
fixed particle number N = |A|/2. Eq. I|15l) is derived by 
considering the leading terms in l|16l) in the limit — » oo 
at fixed A. In the thermodynamic limit Eq. (|15|l means 
that, in contrast to Eqs. JJJJ and lfTU)l. the energy per site, 
e oc (g) 1 has a cusp at g — 1/2, and hence (by Legendre 
transform) the half filled state g = 1/2 corresponds to a 
whole interval of values for the chemical potential. 

For non-zero temperature the energy dependence on 
k may show some rounding due to thermal excitations, 
however there is a cusp in the energy per site, when this 
function is viewed on a scale in which |/c|/|A| > e~ r " c , 
with some r < 1. 

Theorem|21is derived using a representation for the ma- 
trix elements of the relevant operators in the basis that 
diagonalizes {S^.}. An important fact here is that the 
matrix elements of e~ liH are positive in this basis. The 
natural expansion for these matrix elements, e.g., via the 
Lie- Trotter formula, yields a functional-integral represen- 
tation in terms of integrals over "space-time" configura- 
tions, which are represented below by u>. The resulting 
functional integral is not only positive, but also reflection 
positive, and we make use of that fact. However, reflec- 
tion positivity is not essential for the qualitative picture, 
as even without it we would obtain similar results with 
only slightly weaker bounds - with / replaced by 0. 

The measure space over which the integration takes 
place is a Cartesian product of the set of initial spin con- 
figurations times the space of configurations of "rungs" , 
linking pairs of neighboring "space time" sites. A rung 
is parametrized by a pair {x, y} of neighboring lattice 
sites and t £ [0, /?]. For the matrix elements of e~^ H be- 
tween states which correspond to a pair of specified spin 
configurations, one naturally finds an integral over con- 
figurations of arbitrary number of rungs, over which we 
integrate with an "ideal gas" -like measure, in which n- 
tuples are summed and integrated over with the weights 



|-y dt\ ■ ■ ■ dt n . Each rung represents a transformation of 
the spin configuration affected by a specific term in the 
Hamiltonian, and the fugacity-like parameter z is the cor- 
responding amplitude, which in the case of the Hamilto- 
nian considered here is z = 1/2. 

It is particularly convenient to express the spin, or 
particle, configuration in terms of the time-lines of the 
'quasi-particles' which are defined through the occupa- 
tion numbers n x = \ + (—1) X S 3 . There are no quasi- 
particles in the configuration that minimizes the poten- 
tial energy, i.e, if there are ^|A| particles that sit on the 
B-sublattice. The presence of a quasi-particle means ei- 
ther the presence of a particle (S 3 = +1/2), if the site is 
even (A-sublattice), or the absence of one (S* 3 = —1/2), 
if the site is odd (B-sublattice). It is easy to check that 
in this representation the operators a£. and a x act as in- 
sertion of a source and, correspondingly, sink of excess 
spin (relative to the potential- minimizing configuration) , 
although the direction in which the excess spin propa- 
gates changes with the parity of x. Namely, creation of a 
particle on an A site (or annihilation on a B site) results 
in a quasiparticle running 'upward in time', whereas a 
quasiparticle running 'downward in time' originates from 
annihilation of a particle on an A site (or creation on a 
B site). 

Proceeding al ong the above lines, as explained in 
greater detail in [l9||. one obtains 



TrV k e~ 0H = / v 1/2 {dio)e- x ^I[v{u) = k] 



(17) 



where to represents a configuration of a family of disjoint 
oriented loops in A x [0, 0], defined with periodic bound- 
ary conditions in 'time' ([0,/3]), whose orientation alter- 
nates with x, being 'up' along A sites and 'down' along 
B sites. For each configuration, |cj| is the total 'vertical' 
length of the time lines in uj, and v(uj) is the total wind- 
ing number in the periodic 'time' direction. The indicator 
function I[v(co) — k] is 1 if the loop configuration u> has 
total winding number k, and otherwise. The winding 
number can also be computed by adding the spin ori- 
entations of the sites occupied by quasi-particles, along 
any 'constant time' cut through the diagram. The mea- 
sure v z {du>) corresponds to integration, with weights zdt, 
over the times at which the jumps to neighboring lattice 
sites occur, and summation over the possible numbers of 
such jumps. In effect, as mentioned above, the integral 
is over an 'ideal gas'-like distribution of the horizontal 
rungs in the diagram depicted in Figure [3 with the fu- 
gacity parameter taking here the value z = 1/2. That 
value is dictated by the Hamiltonian, where one finds 
1/2 in front of the 'hopping term' S+S~. For later use, 
we find it convenient to consider the measures v z (duj) for 
general z > 0, not only z — 1/2. 
Likewise, for x^y, 



Traia y V k e-P H 



duJ — S(- 



(18) 



v 1/2 (di 



uj e 



I[v{uj) = k] 
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FIG. 2: Loop gas describing paths of quasi-particles for parti- 
cle number N = \A\/2 — 1. A line on an A site means presence 
of a particle, while on a B site it means absence. The hori- 
zontal rungs correspond to hopping of a particle. 



where dui is the set of sources of u>. More precisely, the 
configurations that contribute to the last integral have 
exactly one curve starting at x and ending at y, both at 
time 0, for which we shall use the symbol 7, and other- 
wise only closed loops, of the kind which appear in the 
corresponding trace without the source/sink operators. 

Denoting by A° the set of configurations uj that con- 
sist only of loops, i.e., closed curves, and by A ( ^ y) the 
set of configurations containing one curve connecting x 
with y (at time 0) and otherwise only closed curves, the 
above representation yields for the thermal average un- 
constrained by the particle number 



_ /-AC" 



{sis-) = 



J A o v 1/2 (duj)e-^\ 



(19) 



Now let B^' y) C A { *- y) denote the set of w's which 
consist of exactly one curve connecting x and y, no 
other curves. For a collection of disjoint curves {7^}, 
let A°i 1 Cv4° denote the set of uj's that avoid the col- 
lection {7;} or, in other words, that are consistent with 
it in the sense that insertion of the curves into uj would 
still give an admissible configuration of non-intersecting 
curves. The measures v z (duj) obviously have the product 
structure 



.A(*,y> 



v^e-^ / v z (dw)e 



where we identified uj with 7 in the first integral on the 
right. The following is a convenient bound on the last 
factor. 

Lemma 3 (Contour bound). For any family of dis- 
joint curves {7, }, 



/ v z (duj) 
°~fi } 



e -AM < e J2 s hi\f I Vz (duj)e- X ^ 1 

A" 



where f is the free energy per site. Furthermore, a simi- 
lar bound holds for the integrals further restricted by the 
condition I[v(uj) — 0] and f replaced by fo (which equals 
-(^|A|)- 1 lnTr:Poe-' 3H j. 

The proof, which is presented in Appendix iBl uses the 
chessboard inequality, which is a consequence of the re- 
flection positivity of the functional integral. As was men- 
tioned already, our main qualitative conclusions do not 
require this result, and would follow already from the 
trivial bound in which / and /o, which satisfy / < /o < 0, 
are replaced by 0. 

Lemma |21 implies that 



(S+S-) < f v 1/2 (drf)t 



-ItKW) 



(21) 



To evaluate such expressions, it is useful to consider the 
quantity 

X (z,A) = sup/ u,(d 7 )e- A M 

x,y Jg('>») 

where x = (x,t x ) and y = (y,t y ) are arbitrary points in 
Ax [0,0. 

t 3 „-\t. 



Lemma 4. Iflzd f£ e~ xt dt < 1, th 
X(z,\)< 



1 - 2zd /(f e- xt dt ' 
Furthermore, for any {x, y} and ^, a > 0, 

,(d 7 )e- A ^l <e-«' x - y lx(ze«,A), 



(22) 



(23) 



s (d 7 )e- A|7l /[|7l >*] <e- at x(z,X-a). 



(20) 



Proof. The first inequality here is a random walk bound, 
which is derived by the following 'renewal-type'argument: 
split the integral into a part that comes from curves that 
do not jump at all (which only occurs if x = y) and a 
part where 7 has at least one jump. The contribution 
from the path that does not jump is at most 1. The first 
jump can be in 2d possible directions, hence one gets 

rP 

X<l + 2zd / e" xt dt\. 
Jo 

Since in a finite volume x < 00 a priori, this yields 122fl . 

The rest follows by fairly direct arguments, noting that 
v ze z(d , y) = v z {d^)e^^\ where #(7) denotes the num- 
ber of jumps of 7, which is greater or equal to |x — y| 
in the case in question. The second inequality in ilL'ol) 
is obtained by estimating /[I7I > t] < e ^ 7 ' - *) for any 
positive a. □ 

We shall now use the above functional representation 
to derive Theorem |5J 
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By applying the bounds (|21|l , ill'L'l) and (|23|) to the ran- Defining 
dom walk representation (|19f) . we see that under the con- 
dition stated in (|13[1 the two-point function decays expo- 
nentially, as claimed in l|14|) . 

To prove we start by noting that in case N ^ |A|/2 
the integral in H17fl is over configurations with a non- 
trivial winding number, v{oS) = k. Each such config- 
uration includes a collection of 'non-contractible' loops 
w ith non-zero winding numbers, 1/j = v{^j) i= 0. 
The total length of the set {7^} is at least (3\k\. We shall 
bound the relative weight of such configurations by using 
the second bound in (|23|) and combining it with an argu- 
ment whose purpose is to control the 'entropy' of such a 
collection of long loops. 

Each non-contractible loop jj can be labeled by a 
starting point x_,- £ A where 7^ crosses the time line, 
and a winding number Vj. We shall actually overcount 
by summing over all possible xy's as starting points for 
the loops (with the only restriction that Xj ^ x, for 
j i), and over all possible winding numbers Uj with 

£>,l>|fc|. 



For a given collection {7j} of non-contractible loops, 
we can bound the integral over the remaining loops by 
an integral over loops with zero total winding number 
that avoid the ji's. Hence, starting from l(T7Jl. we get the 
bound 



TrP fc e 



—0H 



< 



E 



n 



{xi,x 2 ,... }CA 

— X 



v x/2 {d lj )e-^I[v{ lj )^Q] 



xl[EMv)\ ^ |fc|]/ni[{7i,72,...}] 
J A o v 1/2 {dLo)e- x ^I[v{u J )=Q] 

j Ao v 1/2 {duj)e-^\l[v{u) = 0] 

Here / n j denotes the indicator function for having only 
non-intersecting loops. Using the chessboard bound of 
Lemma |21 the last fraction can be bounded above by 
Y[. e f°^' x ^\ with f (/3,X) < the free energy per site 
at N — |A|/2. Applying the bounds of Lemma 0] to the 
integral over n loops ji with given absolute value of the 
winding number, wii = \vi\ , we have, for any a > 0, 



H / , v 1/2 (d 7j )e-WM)hi\ lMlj) \ = m .] 

< x{a) n e~ aPY -i m3 ■ (24) 

Here x{a) — x(l/2, A + \fo((3, A)| — a), which is finite if 
a is not too large. 

To complete the bound, we have to sum the right side 
of (|24|l over all the possible choices of the collection of 
the starting points of the winding loops {xi , . . . , x„}, and 
over all possible winding numbers Vj with \v,j\ > 1 and 
J2j Wj \ > \k\- To do so, we employ the following device. 



|A|/2 

with 8 = e~ a @ , we see that the sum in question is given 
by the sum of all the coefficients of z l in P(z)' A l with 
powers I > \k\. Hence 



TrV k e-P H 1 
TrP e-P H ~ 2^i 



dz 



-1 1 



v J- 



where R can be any number greater than 1. The contour 
integral serves as a filter, selecting for us the relevant 
coefficients of the polynomial P(z)' A '. A simple bound 
shows that the above quantity is bounded from above by 



inf 



1 



K-R<l/<5 1 — R 



3 |A|x(«)T^fe 



We now choose R = k / (5x(a)\A\) , assuming that this 
quantity is greater or equal to 2. Note that 8R < i, 
since < ||A| and x(a) > 1. Hence we obtain 



TrToe-f 311 



< 2e 



-aP\k\ 



e 2 x(«) 



\k\ 



1 \k\ 



This inequality is also valid, however, if |/c|/(5x( a )|A|) < 
2, since the resulting bound then exceeds 1, which is 
greater than the left hand side (as shown in Appendix Ej . 
This proves the claim made in , which presents suffi- 
cient conditions for the existence of a cusp in the energy 
dependence on N, i.e., of a gap in the chemical potential. 



VI. NON-CONSTANCY OF THE DENSITY 

In Section III II above we have demonstrated the exis- 
tence of BEC for small A and T, and also that the con- 
densate wave function is constant. Despite this fact the 
particle density has the periodicity of the external poten- 
tial and is not constant for A ^ 0. More precisely, the 
following result is proved below. 

THEOREM 3 (Non- constancy of the density). Let 

g(x) = 7(x, x) denote the particle density in the ther- 
mal equilibrium state at inverse temperature f3. With 
e(X,0) = |A| _1 (i?) equal to the energy per site, 



1 

W\ 



xeA 



> 



A|e(0,/3)| 2 
2d 2 (3d+ A) ' 



(25) 



We first prove this for the ground state. We will show 
that the ground state energy of H, denoted by E(X), 
satisfies 



£(A)<£(0) + iA|A|-cA 2 |A|, 



(26) 
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with c = iegd~ 2 /(3d + A). (Here e = E(0)/\A\ denotes 
the ground state energy per site at A = 0.) Eq. 
implies (|25[1 by the following argument. Write H = Ho + 
XW, with the obvious notation for Ho and W. Since 
E(X) is a concave function of A, we have 

{W) = £'(A) < g(A) " g(0) < ||A| - cA|A|. 
A 

On the other hand, 
(W) = £ [3 + (-!) X ( 5 x)] = 5 |A| + £(-l) x £(x). 

x£A xGA 

Combining the last two equations, we obtain 



1 

W\ 



£(-l) x £>(x) 



xeA 



> cA, 



which proves our claim. 

It remains to show l|26[) . To do this, let the operator 
C be given by 

C = J (^1) X (S'x'Sy — S^Sy) ■ 

<xy) 

As in the Hamiltonian, the sum is over all nearest neigh- 
bor pairs, each pair counted only once. This operator 
has the nice property that 



VII. THE NON-INTERACTING GAS 

The interparticle interaction is essential for the exis- 
tence of a Mott insulator phase for large A. In case of 
absence of the hard-core interaction, there is BEC for 
any density and any A at low enough temperature (for 
d > 3). To see this, we have to calculate the spectrum 
of the one-particle Hamiltonian — ^A + V(x), where A 
denotes the discrete Laplacian and V(x) = A(— l) x . The 
spectrum can be easily obtained by noting that V an- 
ticommutes with the off-diagonal part of the Laplacian, 
i.e., {V, A + 2d} = 0. Hence 

(-±A -d+ F(x)) 2 = (-±A - df + A 2 , 
so the spectrum is given by 



d±y(Y,i C0S Pi) + A2 ' 



where p e A*. In particular, E(p) - E(0) ~ \d(d 2 



A 2 ) 



-1/2 



p| for small |p|, and hence there is BEC for 



low enough temperature. Note that the condensate wave 
function is of course not constant in this case, but rather 
given by the eigenfunction corresponding to the lowest 
eigenvalue of — ^A + A(— l) x . 



VIII. CONCLUSION 



[C,W] = iH . 

Proceeding as in Ijll jl — l|12 |l . a simple bound of the relevant 
double commutators gives 

e ieC We' ieC <W-eH + £ 2 y |A| 

as well as 

e* c H Q e- isC < H + ie[C, H ] + e 2 ^-|A|. 

Note that the e 2 -terms are of order of the volume, due to 
the fact that both C and H contain only nearest neighbor 
terms. With |0) the ground state of H , we therefore 
have, using (0|W|0> = ||A| and (0|[C,ff Q ]|0) = 0, 

E(X) < (0\e isC He- ieC \0) 

< E(Q)(1 - eX) + ±A|A| + ±e 2 d 2 {3d + A)|A|. 

Now the optimal choice of e is e = 2cX/eo, which finishes 
the proof of 

A similar argument works at positive temperature. It 
shows that the free energy depends non-trivially on A, 
and by the same concavity argument as above this im- 
plies the non-constancy of the density, also at positive 
temperature, as claimed in l|25|l . 



We have introduced a lattice model, which is similar 
to the usual Bose-Hubbard model and which describes 
the transition between Bose-Einstein condensation and 
a Mott insulator state as the strength A of the optical 
lattice potential is increased. While the model is not sol- 
uble in the usual sense, we can prove rigorously all the es- 
sential features that are observed experimentally. These 
include the existence of BEC for small A and its sup- 
pression for large A, which is a localization phenomenon 
depending heavily on the fact that the Bose particles in- 
teract with each other. In the Mott insulator regime we 
prove the existence of a gap in the chemical potential, 
which does not exist in the BEC phase and for which the 
interaction is also essential. Bounds on the critical A as 
a function of temperature are included. 
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APPENDIX A: HALF-FILLING AND 
REFLECTION POSITIVITY 

In this appendix we will show that H has a unique 
ground state which has particle number |A|/2. For A = 0, 
this was previously shown in [2(j. We also establish 
the corresponding result at positive temperature, namely 
that the canonical partition function is maximal for par- 
ticle number N — |A|/2, although we do not prove that 
the maximum is obtained only at half-filling. 

The operator H commutes with ^ x «S X . By a Perron- 
Frobenius argument the ground state of H restricted to 
the subspace with fixed value of X^x^x * s unique. We 
claim that the absolute ground state of H corresponds 
to the value ^] x Sj = 0. To prove this we shall use 
reflection positivity. 

We divide the lattice in a left part and a right part 
A = Al U Aji of equal size. We shall identify the space 
Hi = <S> xG a l C 2 with the space (8> x eA n by identify- 
ing factors reflected in the middle plane. We may there- 
fore write the total Hilbert space as H = Hi ® Hi. We 
may then write 

H = H L ®I + I®H R -\ (S+S y + S X S+), 

(xy)eM 

where Hz and Hzi act on Hi and M denotes the set of 
bonds going from the left sublattice to the right sublattice 
(note that because of the periodic boundary condition 
these include the bonds that connect the right boundary 
with the left boundary). Note that Hz ^ Hr. 

We now change H to the unitarily equivalent operator 
H' for which at all sites on the right sublattice we change 
St -> and Si -> -Si. We have 

H' = H L ®I + I®H L -\ Y, (SiS+ + S x S y ). 

(xy}eM 

The same unitary will change ^ x to S' — S®I—I®S, 
where S = J2 x eA L S* acts 011 ^i- 

Let \tp) € Hi ®H\ be a normalized absolute ground 
state for H 1 with S'\ip) = m\ip). We want to show that 
m = 0. We may write \tp) = ^ n \ipn)> where S®I\ip n ) — 
n\ip n ) and I ® S\ip n ) = (n — m)\ip n ). 

We may consider any state \(f>) <^H\®Hi as an oper- 
ator from the Hilbert space Hi to itself, in the following 
way. We introduce a basis \X) in Hi indexed by subsets 
X C Az , labeling the state with all spins up at the sites in 
X and down elsewhere. We shall refer to this representa- 
tion as the standard basis. Then \<f>) may be represented 
as a function associating a complex number (j)(X, Y) to 
any pair of subsets X, Y C Az, namely 4>(X, Y) is given 
by the inner product of |X)®|Y~) with \(f>). Hence \4>) may 
be identified with the operator <fi defined by the matrix 
elements (X\<j>\Y) = <f>(X,Y). 

If A is an operator on Hi then A ® I\ip) — Aip and 

I g) A\ip) = %/jA t where A T is the transposed operator 
represented by the matrix A T (X, Y) = A(Y,X) in the 



standard basis. (Note that transposition is not a canon- 
ical operation, but depends on the basis in which it is 
defined.) 

The operator S is represented by a real symmetric 
matrix in the standard basis. Thus in the above rep- 
resentation ip n maps the subspace where S = (rz. — m) to 
the subspace where S — n and vanishes on the orthog- 
onal complement. Hence we see that tp^tp = J^n^n^n 
and likewise ipip< = J2 n ^ n ^n- ^ follows from this that 
Sft$ = fo$S and £$0t = $$S. Hence 

5(^)1/2 = (^t^)l/2^ s{^) 1/2 = {^ft) 1/2 S. 

(Al) 

Let \tpi), ^2) G Hi <8> Hi denote the states such that 

i>l = (^t)l/2 and ^ = (^t^)l/2. Then |^ x ) and |^ 2 ) 

are normalized (since (V'iIV'i) = Tr (ip\vpi) = Tr (tpft) — 
(ip\ip)) and IHJ implies that S'tyi) = and S'\ip 2 ) = 0. 
We shall prove that 

{«I>\H'W) > \&i\H'Wi) + \(ik\H'\ih)- (A2) 

Since \ip) is an absolute ground state we see that \ipi) 
and ^2) are also absolute ground states. Since they both 
have S' = and the ground state with this property is 
unique we conclude that \ipi) — \1jj2), i-e. , ip^ip — ^ip* . 
Then since we are using a representation in which the 
matrix for 5* is real and symmetric we have 

m = (ip\S'\i/}} = (i/>\(S®I-I®S)\'ip) 
= Tr $$S) - Tr ($$S) = 0. 

It remains to show the reflection positivity (|A2() . We 
may rewrite 

(ip\H'\ip) = Tr $$H L ) + Tr (ft$H L ) 
~ \ E ( Tr + Tr (ftS-$S+j) . (A3) 

Here Mz denotes the set of sites in A^ that connect to 
a bond in M, i.e., sites in A^ that are nearest neighbor 
to a site in A_r. We have used that the operators are 
represented by real matrices in the standard basis. 
The inequality l|A2f) now follows from the inequality 



Tr (ft ArjiA*) < [Tr ($ft) x l 2 A$ft) x l 2 A^ 

r {[ft^f^A^i)) 1 / 2 ^ 



1/2 



1/2 



which holds for any operator A. This inequality is a sim- 
ple application of the Cauchy-Schwarz inequality if one 
uses polar decomposition, i.e., the existence of a partial 
isometry U such that 

xp = U{ft^) 1/2 and {i)ft) 1/2 = Uifti)) 1 ' 2 ^ . 
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At positive temperature we may consider the parti- 
tion function for H restricted to the subspaces with fixed 
value of J2 X We define 

Z(m)=TcV m exp(-0H), 

where V m is the projection onto the eigenspace of J2 X <$x 
corresponding to the eigenvalue m. We claim that the 
partition function is maximal at half-filling, i.e., 



Z(m) < Z(0). 



(A4) 



To prove this we shall again use reflection positivity. 

We first note that the unitary change which mapped 
H to H' will take V m into the operator 



where P n is the projection operator in TL\ projecting onto 
the eigenspace of S with eigenvalue n. Observe that P m 
is a real matrix in the standard basis. 

Using the Trotter product formula we may now write 



Z(m) = lim TrV^x 

k — >oo 

If we use that 



3 2i^(xj)£M D x J y 



xy ■ :< 



= II 



we see that the trace above may be written as sums of 
terms of the form A n A n - m , where 



A n = Tr ni ( P n e-^l k T x e-^l k T 2 ■ ■ ■ e^ H -l k T, 



and each of the operators T\, T2 . . . is a monomial in the 
variables (/3/2fc) 1 / 2 S*±, x 6 M L . 

Since A n is real for all n we see that A n A n ^ m < A n /2+ 
A n _ m /2. If we insert this above and simply undo the 
calculation we arrive at l|A4() . 



APPENDIX B: REFLECTION POSITIVITY 
CONTOUR BOUND 

In this appendix we derive Lemma [3 using reflection 
positivity arguments. 

The measures v z {du)) are reflection positive in the 
following sense. Draw a hyperplane, cither vertically 
(through bonds) or horizontally, that divides I = A x 
[0,0] into two congruent parts U Tr. For any con- 
figuration u>, let uj be its natural reflection through the 
hyperplane (reversing its direction), and for any function 
h on the space of configurations let h(u>) = h(oj). Then, 



for any such complex valued function that depends only 
on the restriction of to to T^, 



h(w)h{w)* e~ XM v z {dw) > 0. 



(Bl) 



This can be seen by noting that once the behavior of co 
on the hyperplane is fixed, the distribution of the left 
and right sides (or top and bottom) are conditionally 
independent, and are mirror images of each other. 

Reflection positivity leads to what is known as the 
chessboard inequality |2lll22|. In essence, it is a multiply 
reflected generalization of the Schwarz inequality, which 
allows us to obtain bounds on the expectation value of 
a product of local variables in terms of thermodynamic 
quantities. 

The function whose average we need to estimate is 
X£>(k>) - the indicator function which is 1 if the curves in 
u> avoid a specified set Del and otherwise. One may 
start by partitioning the imaginary time interval [0, /3], 
and correspondingly the "space-time" T, into equal short 
segments whose reflections tile T. For any subset D C T 
that is a union of elements of the finite partition of Y the 
strategy, which is explained in detail in [2l]], yields 

J Ao v z (duj)xD(<jj)e- x ^ 
J Ao v z (duj)e- x ^\ 



< 



J Aa v z (duj)xr(uj) 
£40 v z (dto)e- x 




D\/\T\ 



e'W, (B2) 



where \D\\s the total length of D. By refining the parti- 
tion, and applying elementary continuity arguments (the 
dominated convergence theorem) , we conclude that l|B2|) 
extends to all sets D which are finite unions of closed 
intervals. This proves 12UI) . 

To prove the second statement in Lemma |21 we note 
that reflection positivity holds also for the restriction of 
the measure to lo's with winding number. I.e., for h as 
in ifEST) . 



h{ui)h{uj)* I[v{uj) = 0]e- x ^v z {dLj) > 



If the hyperplane dividing T into UTj; is horizontal, 
this is clear, since fixing w on this hyperplane fixes the 
winding number. If it is vertical, however, we note that 
I[v{ui) = 0] can be written as J^ fe Ik(uj)Ik(u>), where k 
runs from to |A|/2, and Ik is the indicator function 
for to restricted to Tl having winding number k. (Note 
that to restricted to may have sources and sinks on 
the boundary, and when counting the winding number 
we also have to consider the curves resulting from these.) 
The second claim made in Lemma[!|]follows by proceeding 
as in i|B2|) . but with the added restriction to winding 
number. 
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